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We explore a strategy for protecting the evolution of a qubit against the effects of environmental 
noise based on the application of controlled time-dependent perturbations. In the case of a purely 
decohering coupling, an explicit sequence of control operations is designed, able to average out the 
decoherence of the qubit with high efficiency. We argue that, in principle, the effects of arbitrary 
qubit-environment interactions can be removed through suitable decoupling perturbations acting on 
the system dynamics over time scales comparable to the correlation time of the environment. 

I. INTRODUCTION 

Decoherence remains one of the most serious obstacles to the exploitation of the speed-up promised by quantum 
computation jD. Broadly speaking, two different philosophies are being investigated to overcome the decoherence 
problem. On one hand, passive error-prevention schemes have been proposed, based on the idea of encoding logical 
quantum bits (qubits) within subspaces which do not decohere owing to symmetry properties |^,||. On the other hand, 
active error-correction approaches have been formalized within a sophisticated theory of quantum error-correcting 
codes (QECC), where a logical qubit is encoded in the larger Hilbert space of several physical qubits and suitable 
feedback operations are conditionally carried out |||| . 

Although a purposeful manipulation is implied in the latter case, quantum error-correcting codes can be properly 
interpreted in terms of a clever redundancy in the software architecture rather than a physical way to operate on 
decoherence. In this work, we explore the possibility of using control techniques to modify and eliminate decoherence. 
Unlike recent proposals for feedback (or closed-loop) control schemes of decoherence in quantum optical systems || , we 
apply control in the simpler open-loop configuration to general models of quantum information processing systems . 
The underlying idea is suggested by high-resolution pulsed Nuclear Magnetic Resonance (NMR) , where astonishingly 
versatile refocusing and decoupling techniques are nowadays available to remove the effects of interactions among 
the spins that are considered unwanted or uninteresting j^]. In our analysis, we outline the conditions under which 
analogous procedures can be extended from eliminating interactions internal to the system to suppressing interactions 
of the system with an external quantized environment. In particular, the role of the environment correlation time will 
be pointed out as a further parameter to be engineered in the struggle for preserving quantum coherence. 

II. QUANTUM BANG-BANG CONTROL OF QUBIT DECOHERENCE 

We start by investigating a prototype situation that conveys the basic idea in the simplest form. We will focus on 
the dynamics of a single memory cell of quantum information (qubit) undergoing decoherence due to the coupling to 
a thermal reservoir. The physical qubit can be associated either to a fictitious or to a real spin- 1/2 system, the latter 
case allowing for a direct reference to the language of Nuclear Magnetic Resonance and NMR quantum computation 
[^|,^) . We assume that the fastest relaxation process originated by the interaction with the quantized environment is a 
purely dephasing process i.e., in NMR terminology, no dissipative Xi-type of decay takes place. A general model for 
the dynamics of the overall qubit + bath system is then provided by the purely decohering spin-boson Hamiltonian 

(h = iy. 



H = H s +H B + H SB = ^-°z+J2 Wfcfe fe fefe + ^ kb l + 9kh) ■ (1) 

k 



Here, o z is the standard diagonal Pauli matrix, with qubit basis states denoted as \i), i = 0, 1, while b\,bk,gk are 
bosonic operators and coupling parameters for the A:-th field mode respectively. Hamiltonian (|l|) is widely used in 
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the quantum computation literature to investigate the effect of phase errors, representing the most nonclassical and 
dangerous source of errors for quantum qubits Since [a z ,Ho] = 0, spin populations are not affected by 

time evolution and decoherence dynamics is characterized completely by the qubit coherence with respect to the 
computational basis: 

PQ1 (t) = <0| p s (t) |1) = (0| Ti B {ptot(t)} |1) = Tr B , s {ptot(to)<T+(t)} = WW) , (2) 

where, starting from the left and using standard notations, the relevant reduced density matrix element poi(t) — Pio(t) 
in the Schrodinger picture is linked to the expectation value of the time-evolved ladder operator in the Heisenberg 
representation 



<>■+(*) =oM*) +**»(*))■ (3) 



The evolution of the qubit coherence (Q) can be calculated exactly under the customary assumptions about the 
initial state of the overall system, i.e. qubit and environment are initially uncorrelated and the environment is in 
thermal equilibrium at a temperature T. In the Heisenberg representation, one may write formally 

(a+(t)) = (G tot (t ,t)a + (t )) , (4) 

with the propagator G tot (to > t) determined by the solution of the Heisenberg equations for the coupled spin + bath 
motion: 



o+(t) = iu; c7 + (t) + 2ij: k (g k bl(t)+g* k b k (t))a + (t) , 
b k (t) = -iw k b k (t) - ig k a z (t) , 
{ b{(t) = +iu k b{(t)+ig* k a z (t) . 



(5) 



The result can be written in the following form j7j]: 

P01 (t) = e ^°(*- t °)- r <> (*"*<>) poi(*o) , (6) 
the loss of phase information being characterized by the damping function 

r (t -to) = £ r (fc; t - 1 ) = £ A\ 9k \> coth fe) , (7) 

k k \ / U k 

in units where the Boltzmann constant k B = 1- In the limit of a truly macroscopic environment, a description in terms 
of a continuum of modes is appropriate and the dependence of the decoherence function (Q) on reservoir properties 
can be cast in a compact form after introducing the spectral density function I (to), 

( UJ \ 1 — cosw(< — to) 



ro(t-t )=J (LjT (cj;t-t )= J ^/M4cothl-] ^. (8) 

Depending on the temperature T and the spectral density I(uj), qualitatively different open-system evolutions arise in 
general, with a different interplay between quantum fluctuation and dissipation phenomena. Regardless the details of 
the spectral density function, however, the existence of a certain ultraviolet cut-off frequency cu c is always demanded 
on physical grounds, leading to 

I(uj) ^0 for uj > uj c . (9) 

Although the specific value of lo c depends on a natural cut-off frequency varying from system to system, u c can 
be generally associated to a characteristic time r c ~ oj^ 1 setting the fastest (finite !) time scale of the irreversible 
dynamics. r c is known as the correlation time of the environment. The dynamics of the decoherence process arising 
from (||) for various choices of has been investigated in detail elsewhere ||],[nj ■ A pictorial representation for the 
important class of Ohmic reservoirs, I(uj) oc we - "/" , is shown in Fig. 1. 

We introduce now a procedure aimed at improving the coherence properties of the qubit by adding a controllable 
time-dependent interaction to the original Hamiltonian: 

H(t) = H + H 1 (t) = H S + H B + H SB + Hi (t) . (10) 
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In the same spirit underlying multiple-pulse techniques in the manipulation of nuclear spin Hamiltonians ||, we 
try to average out the unwanted effects of the qubit-reservoir coupling Hsb by applying a sequence of coherent n- 
pulses that repetitively flip the state of the system. Under the assumptions that the duration of the pulses is short 
enough compared to the typical decoherence time and the strength of the control field is sufficient to override the 
Hsb coupling, the pulsed-mode operation allows one to separate the actions of the bath and the external field, by 
neglecting Hsb while Hi (t) is on. Specifically, H\ (t) is assumed to schematize a train of np identical 7r-pulses along 
the x-axis applied on resonance at instants t — tp"\ n — 1, . . . , rip, with pulse separation tp^ +1 ^ — t^ = At. By 
invoking, as usual, the rotating-wave approximation, we have 



with envelope function 



Hi(t) — S(i) [ cos{u>Qt)a x + sm{u>ot)a 



v J ■ 



(11) 



Z(t) = Y,v[9(t-tP)-#{t-tP-Tp) 



An) 



to + nAt, n = 1, 



,np . 



(12) 



In Eq. (|12|), $(■) denotes the Heaviside step function, and the height V and the width rp of each pulse satisfy 
2Wp = 7r. To simplify things, we work henceforth in the limit of infinitely narrow pulses rp — > 0, assuming the kicks 
of radiofrequency control field large enough to produce instantaneous spin rotations. By analogy with the classical 
technique of bang-bang (or on-off) controls, whereby piecewise controls with extremal values are exploited one 
may look at this strategy as an implementation of quantum bang-bang control. 

In order to depict the evolution associated to a given pulse sequence, it is convenient to think the latter as formed 
by repeated elementary cycles of spin-flips, a complete cycle being able to return the spin back to the starting 
configuration. For definiteness, let us analyze the first cycle, made of the following steps: evolution under Hq during 



t < t < t { p ] 



Tr-pulse P\ at time t^p; evolution under H during tp^ 



<t< tp\ 7r-pulse P2 at time tp\ After a total 



time t\ = to + 2At, the first cycle is complete. The description of 7r-pulses turns out to be extremely simple in the 
Heisenberg representation. Nothing happens to the bath operators bk,b\ in the limit of instantaneous pulses, while, 

by denoting with t p the times immediately before (after) a pulse respectively, spin operators are transformed as 
follows: 



a+(t+) 



Vz(tp) , 
r + (tp)V 



(13) 



In terms of the free propagator Gtot(ti,tj) introduced in (|J) to evolve coherence from ti to tj, the time development 
during the cycle can be represented as 



<o+(*o + 2At)> = (Gtot(t , to + At) a + (t ) G\ ot (to + At, t + 2 At)) , 



(14) 



to be compared with 



<<7+(t + 2At)} = {G tot {to + At, to + 2At) G tot (t , h + At) a+(t )) 



(15) 



in the absence of pulses. Since instantaneous rotations introduce discontinuous changes in operators (13), care must 
be taken in evaluating the two propagators G to t(to,to + At), G to t(to + At, t + 2 At) separately, by solving Heisenberg 
equations of motion with initial conditions at t = to, t = tp~ = to + At respectively. Only at the end of the 
calculation everything can be expressed with respect to the initial time of the cycle. Omitting the details, the result 
for the coherence evolution over the first complete cycle is 



Poi(t Q + 2At) = e~ r ^ N=1 ' At ) p 01 (t Q ) 



(16) 



where a new decoherence function for N = 1 spin cycles has been introduced: 

,l' (fc;2A£) 1-2 

k 



r P {N = l,At) =^r (fc;2At) 



(17) 



Since, for each mode, the additional factor arising from the pulses is of order 0(id\At 2 ) ^ 1 for small At, we may 
guess that something interesting is happening in a regime where the state of the qubit is tipped very rapidly. This 
is made clear by generalizing the description to an arbitrary number N of spin-flip cycles, involving a total number 



3 



of rip = 2N 7r-pulses. After straight forward calculations along the same line outlined above, the expression for the 
qubit coherence at the final time tjy = to + 2N At is the following: 

p 01 (t + 27V At) = e-^W^ p 01 (t ) , (18) 

with 



T P (N, At) = J2 r o(fc; 2NAt) 1 - f k (N, At) , (19) 



fc 



N 



MN, At) = 2 _ * a 8 *- 1 *** . (20) 



n=l 



From inspection of Eq. ( |l9| ) , the contribution due to the pulse sequence turns out to manifest in the typical form of an 
interference factor. The implications for the decoherence properties are easily stated by considering the mathematical 
limit where At — ► 0, N — > oo, subjected to the constraint 2N At = t^ — to. Under these conditions, one can prove 
that 

lira f k (N, At) = 1 Vfc lim T P (N, At) = . (21) 

Thus, in the limit of continuous flipping, decoherence is completely washed out for any temperature and any spectral 
density function. 

Obviously, a continuous limit of this kind is scarcely meaningful from a physical point of view. However, this ideal 
situation should be approached if At is made small compared to the fastest characteristic time operating within the 
environmental noise, i.e. the reservoir correlation time r c . Hence, we expect that a sufficient condition in order to 
meet (|2l]) and suppress decoherence is 

At < t c . (22) 

More explicitly, this result implies that, given an arbitrary time i, one can always recover the initial state and the 
coherence of the qubit by making t the end time of a pulse sequence and by adjusting the parameters to satisfy 
t = tN = 2 N At and At -C t c . Then, at time t, a coherence echo is formed. Alternatively, by keeping the qubit flipped 
and restricting the observation to cycle times f/y, N = 1, 2, . . ., the system is found to evolve ideally as it would do in 
the absence of the coupling Hsb responsible for decoherence. A typical behavior originated by the pulsing procedure 
for the prototype high-temperature Ohmic environment of Fig. 1 is displayed in Fig. 2. 

So far, the suppression of decoherence has been derived in a rather formal way. Actually, a simple physical 
explanation can be provided as well. Similarly to the original well-known spin-echo phenomenon Jl4| , and to the more 
sophisticated solid-echoes or magic echoes experiments ||, the basic argument here is a time-reversal argument. The 
examination of a single elementary spin-flip cycle suffices to capture the underlying mechanism. Roughly speaking, 
and looking back at the representation (|J), it is the presence of the transformed propagator G^ t (t + At, to + 2 At), 
generated by the couple of 7r-pulses, that simulates the effect of a time-reversal. Would the evolution during the 
second half of the cycle be identical to the one in the first At interval, then it would be 

G t "t(to + At, t + 2 At) = G to t(to, t + At) , (23) 

and, therefore, (er+(to + 2 At)) = (er + (to)) as a consequence of the cyclic property in the trace. Instead, this reversal 
is only approximate in general since the two propagations differ by a dephasing factor e lWfeAt in the evolution of 
each reservoir mode [Q. However, if the condition (|2^) is met, then the cycle is effectively equivalent to an exact 
time-reversal and, by iteration on every cycle, the elimination of decoherence ( |2l| ) is achieved. 

III. DYNAMICAL DECOUPLING OF QUBIT-ENVIRONMENT INTERACTIONS 

In this section, we rederive the result established above in a form that opens up the way to further generalization. 
The first step is to formally reinterpret the method of reducing environment-induced decoherence by successive 
application of 7r-pulses within the general framework of decoupling techniques based on controlled averaging. In 
NMR, sophisticated decoupling schemes arc routinely used to simplify complex spectra by manipulating the underlying 



spin Hamiltonian to an extent allowing for a successful analysis [p 15 1. In particular, a relevant class of decoupling 



procedures, including spin-decoupling and multiple-pulse experiments, involves selective averaging in the internal spin 
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space. The idea is to introduce controlled motions into the system, with the time-dependence designed in such a 
way that undesired terms in the Hamiltonian are averaged out. In extending similar techniques to the decoupling 
of interactions between a system and its environment, the major difference stems from the fact that the decoupling 
action can be easily exerted only on the system variables, the bath degrees of freedom being generally uncontrollable. 

Looking back at the Hamiltonian we start by seeking a perturbation H\(t) to be added as a suitable decoupling 
interaction in order to remove Hsb, Eq. (|Io|). We restrict to a situation where Hi(t) is cyclic, i.e. satisfying the 
following conditions (to — henceforth): 

(£) Hi{t) =Hi(t+ At) for some At; (24) 
(ii) Ut(t) = Tcxp j- ij dsi?i(s)| = Ut(t + T C ) for some T c . (25) 

From (ii), U\(T C ) =1 and T c is called the cycle time. An elegant description of the dynamics arising in the presence 
of a cyclic perturbation is provided by the so-called average Hamiltonian theory [plJlql . We only recall here the basic 
ingredients in a language that is appropriate to quantum information processing. In particular, we schematize the 
perturbation as a sequence on np ideal pulses (or logic gates), represented by unitary transformations Pi,... ,P np , 
and separated by free evolution periods under H . The operators {Pi} fulfill Ui(T c ) = P np . . . P2P1 = 1 by cyclicity. 
Average Hamiltonian theory is based on moving to the time-dependent interaction representation with respect to 
Hi(t) defined, as usual, by 

Ptot(t) = Ui{t)p tot {t)U\{t) , (26) 

with 

J t Ptot{t) = -i[H{t),p tot {t)\ , H(t) = ul(t)H Q Ui{t) . (27) 

In (|27|), Ho is given by (0) and, according to the standard NMR literature, the transformed Hamiltonian H(t) is also 
known as the toggling frame Hamiltonian. If Ut„±(t) and Utot(t) denote the time evolution operators in the Schrodinger 
and interaction picture respectively, due to ((2^) one gets 

U t ot{nT c ) = Ui{nT c )U tot {nT c ) = U tot {nT c ) = [U tot (T c )\ = e~ lH *" , (28) 

i.e. provided the observation of the dynamics is restricted to stroboscopic and synchronized sampling, t n = nT c , it 
is sufficient to know the evolution over a single cycle as described by the transformed Hamiltonian H(t). In (|2^), 
the last equality defines the average Hamiltonian and contains the main result of average Hamiltonian theory: for 
suitable repeated observations, the motion of the system under the influence of the time-dependent field Hi(t) can 
be represented by a constant average Hamiltonian H. The calculation of H is usually performed on the basis of a 
standard Magnus expansion of the time-ordered exponential defining Utot(T c ) |lq1 , i.e. 

U tot (T c ) =Texp| - l J^ dsH(s)^e-^ T ^e-^ m + 7lll) +-^ . (29) 

In our case, the evaluation of ([2?]) is simplified since the toggling frame Hamiltonian H(t) is piecewise constant during 
the intervals Atk separating consecutive rotations and one can introduce stepwise transformed Hamiltonians: 

Utot(T c ) = e -^"p A *-p . . . e -^° At ° , H k = {P k ... P^ 1 H (P k ... Pi) , (30) 



where, obviously, ^ fe At^ = T c . Thus, the lowest-order approximation H^ to the average Hamiltonian in (29) has 
particularly simple form: 



I f If np 

H () = — H At + ... + H np At np I = — J2 At k Pi 1 ■ ■ ■ Pu 1 H Pk...P 

±c L > Jc fe=o 



(31) 



In addition, explicit expressions for the corrections arising from higher-order terms are systematically available from 
the Magnus expansion |f5| . Loosely speaking, the goal of decoupling is to devise a tranformation to a toggling frame 

(p6|), where the unwanted coupling Hsb no longer appears up to a certain order H r and higher-order contributions 

H , . . . are made neglegible. 
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We are now in a position to apply this formalism to the evolution of the decohering qubit, whose Hamiltonian we 
rewrite for convenience as follows: 

Ho=H 8 +H B +H 8 B = Y (r * + ^2 U)kb k bk + <T * B '> ^ = J2^kbl+9* k b k ) . (32) 

k k 

Using (p7|), the transformation to the toggling frame associated to H\(t) leads to 

H(t) =H B + u\{t) H s U x {t) + [u\{t) a z U^t)} B z , (33) 
and to a zero-th order average Hamiltonian given by (|3l|): 



i np r 1 np 

H (Q) =H B + 7 F J2 ^tkV^HsVk + — ^ ^kV^o z V k 



c k=0 



B z , (34) 



where, in (|33|) and (|3J), the fact that reservoir operators are unaffected by the control field has been evidenced and 
the short notation Vk = Pk ■ ■ ■ Pi has been introduced. The second and third terms in ( p4[ ) correspond, in general, to 
a transformed qubit Hamiltonian and a transformed qubit-bath interaction. It is immediate to realize that the effect 
of the 7r-pulse sequence of Sec. II viewed in this frame is to cause the latter interaction term to vanish. To make the 
identification explicit, we rearrange the elementary spin-flip cycle as follows: 

At 18n At At lsn At At , . 

-T- p - -T-T- p * -T-T- (35) 

cycle 

where the rotation axis of the pulses has been indicated and, to compare with (|3l]), Ato = At2 = At/2, At± = At, 
T c = 2 At. Written in the form (|35|), the decoupling sequence for ( |32| ) is nothing but a variant of the famous Carr- 
Purcell (CP) sequence, that is ordinarily exploited to get rid of static applied-field inhomogeneities |§,|l5|. At variance 
with this standard usage of the CP-sequence, however, where the size of At is of no importance, the averaging of Hsb 
at zeroth-order does not guarantee by itself the elimination of decoherence. Obviously, higher-order terms in (|29| ) 
have to be quenched. The question under which circumstances the Magnus series can be truncated after the leading 
term or the few lowest-order corrections is nontrivial. It is possible to show [|l6| that, since the r-th order contribution 

— 0(At r ), a sufficient condition is Atuj c <§; I, which is identical to (|22|). Then, from (f28|), one ideally gets the 
stroboscopic equality 

Ps (nT c ) = Tr B {e"^ (0) "^ p S (0)p B (0) } = Ps (0) . (36) 

The generalization of the above scheme to decouple arbitrary qubit-environment interactions is, in principle, straight 
forward since the most general bilinear coupling can be expressed as a mixed sum of error-generators |2],|l7| : 

Hsb = cr x B x + (T y By + <j z B z , (37) 

for suitable reservoir operators. Then, provided there is no constraint on the rate of control so that condition ( [22] ) 
can be assumed, one has to ensure the existence of a gate sequence generating the required temporal average: 

np 

J2 At k VkO a Vk = , a £ {x, y, z} . (38) 

fc=0 

The purely decohering coupling corresponds to B x = B y = 0. Actually, binary sequences of the Carr-Purcell type also 
suffice to decouple any form of interaction ( |37j ) involving at most two error generators. A specially relevant case is a 
Jaynes-Cummings-like dissipative coupling with B z = and ^ 

H SB = J2 (9kb\a- + h.c.) => B x = J2 (fffc&l + h.c.) , B y = Hs*4 + h.c.) , (39) 

which can be eliminated, in principle, by a sequence of 7r-pulses along the i-axis. A slightly more elaborated sequence 
is necessary to decouple the qubit from the simultaneous action of the three error generators. It turns out that it can 
be derived on the basis of a simple group-theoretic argument. A more detailed and general formulation of the method 
is presented elsewhere [M . 
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IV. CONCLUSIONS 



Our work demonstrates the possibility to modify the evolution of a quantum open system by applying external 
controllable interactions. From the perspective of quantum information, the analysis suggests a different promising 
direction compared to conventional quantum error-correction techniques. The practical usefulness of the proposed 
approach strongly depends, in its present status, on the time scale of the motional processes causing relaxation. The 
effectiveness of analogous schemes under less idealized assumptions and in the presence of a finite bound on the control 
rate deserves further investigation, together with the possibility of examining decoherence properties within a fully 
quantum mechanical control configuration as recently proposed in |lq| . 
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FIG. 1. Qubit decoherence as a function of time for an Ohmic environment. Time is in units of T _1 and lo c = 100. High- 
and low-temperature behaviors are depicted, (H) cj c /T = 10 -2 and (L) w c /T = 10 2 respectively. 



FIG. 2. Qubit pulsed decoherence as a function of time for the Ohmic high-temperature environ- ment of Fig. 1. A pulse 
separation At = r c /10 has been used and coherence stroboscopically evalua- ted using Eqs. (18)-(19). Long-time deviations 
from unit value arise from cumulation of errors in the presence of a small but finite At. 
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